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Abstract. We predict the survival time of initially bound star clusters in the solar neighbourhood taking into 
account: (1) stellar evolution, (2) tidal stripping, (3) shocking by spiral arms and (4) encounters with giant 
molecular clouds. We find that the predicted dissolution time is tdis = 1.7(Mi/10* Mq)° ®'^ Gyr for clusters in the 
mass range of 10^ < Mi < 10^ Mq . The resulting predicted shape of the logarithmic age distribution agrees very 
well with the empirical one, derived from a complete sample of clusters in the solar neighbourhood within 600 
pc. The required scaling factor implies a star formation rate of 4 10^ Mq Myr~^ within 600 pc from the Sun or a 
surface formation rate of 3.5 10~^° M0yr~^pc~^ for stars in bound clusters with an initial mass in the range of 
10^ to 3 10* Mq. 
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1. Introduction 

' The first empirical determination of t he life t ime o f clus- 
. ters in the solar neighbourhood is by lOord (Hm, who 
' noticed the lack of clusters old er than a few Gyr in 
the solar neighbourhood. Later, IWielerJ l)l97l|) derived 
a mean dissolution time of 0.2 Gyr from the age dis- 
tribution of clusters. Since most of the observed clus- 
ters within about 1 kpc from the sun have a mass in 
the range of lO'^ to a few 10^ Mq the value derived by 
Wielen is for clusters in that mass range. Theory pre- 
dicts that the dissolution time of clusters depends on 
their initial mass in that massive clusters survive longer 
than low mass clusters (e.g. ..Spitzer 1958, Wieleu 198^, 
iGhernoff fc WeinberdligMlGnedin fc Ostrikedll997l and 
references therein). Bau mgardt fc Makinol l)2003^ (here- 
after BM03) showed from TV-body simulations that the 
dissolution time of clusters in the tidal field of the galaxy 
depends on their initial mass. Mi, as tdis ~ M°-^^. 
Independently, this same power-law dependence was also 
derived em pirically in a study of cluste r samples in four 
galaxies bv lBo^itloukos fc LamersI l(200.ll) . 

The dissolution time of clusters in th e solar neigh- 
bourh ood was recently redetermined by iLamers et alJ 
('2005') (hereafter L05), based on a new cluster sample of 
[k ■larciiciikd et al- (JOO^- They found a dissolution time 
of tdis = t4iMi/ 10* Mq)°-^^ with t4 = 1-3 ± 0.5 Gyr for 
clusters with 10^ < M < 10*. This is a factor 5 shorter 
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than the t^ — 6.9 Gyr that follows from the tdis{Mi) re- 
lation derived by the A^-body simulations of BM03 for 
clusters more massive than 4500 Mq at a Galactocentric 
distance of 8.5 kpc. The simulations of BM03 include a 
realistic stellar mass function, stellar evolution, two-body 
relaxation, a detailed treatment of binary evolution and 
close encounters of stars. Part of this difference may be 
due to the fact that in low mass clusters, with lifetimes 
shorter than about 1 Gyr, the dynamical evolution is af- 
fected by the presence of massive stars during most of their 
lifetime (see Fig. 5 of BM03). It is doubtfuU that this effect 
alone can fully explain the difference between the results 
of L05 and BM03. In fact, the large discrepancy suggests 
that other, probably external, disruptive effects must play 
an important role in destroying star clusters in the solar 
neighbourhood. 

In this paper we explain the lifetime of clusters in the 
solar neighbourhood, by taking into account the combined 
effects of stellar evolution, tidal stripping, encounters with 
giant molecular clouds (GMCs) and spiral arm shocks. 
We use stellar population models to describe the stellar 
evolution and the results of BM03 for tidal stripping. For 
the effects of GMCs and spiral arms we adopt the new 
estimates from the recent studies bv.Gieles et al. (,2(3063) 
(hereafter GPZB06) and iGieles et alJ |2006ah (hereafter 
GAPZ06), which arc based on TV-body simulations. 

The structure of the paper is as follows. In Sect. 2 
we discuss the predicted mass loss from star clusters by 
stellar evolution, tidal stripping, encounters with GMCs 
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and spiral arm shocks. We calculate the mass evolution of 
clusters due to these four effects. In Sect. 3 we compare 
our results with the observed age distribution of clusters in 
the solar neighbourhood. The discussion and conclusions 
are given in Sect. 4. 

2. The decreasing mass of star clusters 

2.1. Mass loss by stellar evolution 

The mass loss from clusters due to stellar evolution has 
been calculated for cluster evolution models by several 
groups. We adopt the GALEV models for single stel- 
lar populations with a Salpeter tvp e mass function in 



the range of 0.15 < M/Mp, 
lAnders fc Fritze-v. Alvenslebe: 




kilz et alJ l2002t 

These models are 
based on stellar evolution tracks from the Padova group, 
which include ov ershooting, mass loss due to stellar wind s 
and supernovae. ijBertelli et al ]ll994H rxirardi et alJ l20n(l(l . 
L05 have shown that the fraction of the initial cluster mass 
that is lost by stellar evolution, gov(i) = AAZ/M^, can be 
approximated accurately by 

loggcv(t) = (logi - a^yf'^- + c^v for t > 12.5 Myr. (1) 

with t in yrs. For solar metallicity models with a Salpeter 
mass function Ocv = 7.00, 6ev = 0.255 and Cov = —1.805. 
This function describes the mass loss fraction of the mod- 
els at t > 12.5 Myr with an accuracy of a few percent. The 
mass loss at younger ages is negligible because stars with 
M* > SOM© hardly contribute to the mass of the cluster. 
The mass loss of a cluster by stellar evolution is 



dM 
dt 



= -M{t) 



cvol 



dt 



(2) 



2.2. Mass loss by the galactic tidal field 

BM03 have calculated a grid of A^-body simulations of 
clusters in circular and elliptical orbits in the tidal field 
of a galaxy for different initial cluster masses, galacto- 
centric distances R, and different cluster density profiles. 
The stars follow a Kroupa initial mass function and stel- 
lar evolution is tak en into account during the evolution. 
iGieles et al.l l|2004) have shown that for all models with 
clusters with Mi > 4500 Mq of BM03 the dissolution time 
can be expressed as a function of the initial cluster mass 
as 



0.62 



(3) 



where is a constant that depends on the tidal field 
strength of the galaxy in which the cluster moves and 
on the ellipticity of its orbit. (If the tidal field were the 
only disruptive process then <4 would be the lifetime of 
a cluster with Mi = 10^ Mq.) The mass loss due to the 
Galactic tidal field can then be written as 



fdM\ 



'M{t) -{M/10*Mq) 



0.38 



tidal 



^dis 



t4/104 



MQMyr-i(4) 



We adopt the value of = 6.9 Gyr from BM03 for clus- 
ters in circular orbits at Rq = 8.5 kpc, although this may 
slightly overestimate the disruption time of the lower mass 
clusters by about a factor 2 or so (see §1). 

2.3. Mass loss by spiral arm shocking 

GAPZ06 studied the dissolution of star clusters by spi- 
ral arms by means of A^-body simulatio ns. They used 
and a djusted the analytical expression of lOstriker et alJ 
lll972ll for the dissolution time of star clusters due to disk 
shock, to derive an expression for the dissolution time 
of star clusters by spiral arms (isp)- Mass loss by spi- 
ral arm shocks will occur just at the moment the cluster 
crosses the spiral arm. Assuming that spiral arms move 
with a constant pattern speed (Qp) and that the mat- 
ter in the disk has a constant circular velocity (Vdisk), 
the relative velocity between the two (Vdrift) depends on 
the location in the galaxy (R). Density waves that pass 
with a low velocity hav e a la rge effect on the star clus- 
ters (e.g. Gjtrikcr ct al] ll972l) . Therefore, the disruptive 
effect of spiral arm shocks is most important close to the 
corotation radius (Rcr), i-e. the point where the disk 
and the spiral arms have the same rotational velocity. 
We adopt the "average" spir al arm model of GAPZ 06, 
which is based on the study of lElmeereen et all l)l989|) of 
the spiral galaxies M81 and MlOO, to derive the densitv 
contrast of the spiral arm. IPias fc Leoinel l|2005(l found 
rip = 25.9 km s~^ kpc~^ for the spiral arms in the Galaxy, 
from a study of the nearby star clusters, and a corota- 
tion radius {Rcr) almost coinciding with the solar ra- 
dius Rcr/Rq = 1.06 ± 0.08. Based on the adopted values 
of i?o = 8.5 kpc, Vdisk — 220 km s^^ and the assump - 
tion that our Galaxy has 4 spiral arms lIVallefjEooi . 
GAPZ06 used Vdrift = 12.5 kms'^ = 12.7pcMyr-i and 
idrift = 1.05 Gyr. Taking into account the ratio / = 
{AM/M)/AE/E) = 0.3 between the energy gain and 
the mass loss, predicted by GAPZ06 we find for the solar 
neighbourhood that 



tsp = 20 



M 



104 



3.75 pc 



= 20 {M/W'^Mq 



1-3A 



Gyr 



Gyr, 



(5) 



where we have substituted the observed m ass-radius rela - 
tion of clusters in nearby spiral galaxies of lLarsenI l)2004j) : 
Th = 3.75 (AZ/IO^Mq)^, with A = 0.10 ± 0.03. The mass 
loss of clusters due to spiral arm shocks is then 



/dM 

[~dr 



-M{t) 

tds 



= -0.5 



( Mjt) 
VlO^Mg 



3A 



Mo Myr- 1.(6) 



Notice that for A = 0.1 the mass loss due to shocking 
by spiral arms has almost the same mass dependence, i.e. 
oc M°-^, as the mass loss by the tidal field (viz. M"-^®). 
We will use the mass loss relation Eq. |B1 as a statistical 
mean. 



H.J.G.L.M. Lamers & M. Gieles: The destruction of cluster in the solar neighbourhood 



3 



2.4. Mass loss by giant molecular cloud encounters 

GPZB06 studied the encounters between GMCs and clus- 
ters with A'^-body simulations. They derived an expression 
for the energy gain and the resulting mass loss for the full 
range of encounter distances, from head-on to distant en- 
counters. Adopting a mean GMC density in the galactic 
plane near the sun of p„ = 0.03Mf7> pc~'^, a surface den- 
sity of GMCs S„ = 170Mq pc'^ |solomon et-ailll987|l 
and a mean velocity dispersion of clusters and GMCs of 
~ 10 kms^^, they derived a dissolution time (tcMc) for 
clusters by GMC encounters in the solar neighbourhood 
of 

r.r. f M \ /3.75pc\^^ 
iGMC = 2.0 ( ttttt:^ ) ( — ^ ) Gyr. (7) 



10000 I 



EVOLUTION 



TIDAL FIELD 



6000 - 



a 4000 



2000 - 



SPlMTTTOfflS 



REMAIKIKG MASS 



GMCs 



5.0x10" 



If we assume the same mass-radius relation as before we 
find that the mass loss rate due to encounters with GMCs 
is 

3A 



1.0x10" 
t [yr] 



1.5x10" 



\~dr 



-5.0 



M{t) 
104 Mp 



M0Myr-i.(8) 



GMC ^GMC K^'^'^v^Q, 

Notice that the mass dependence is the same as for dis- 
solution by spiral arm shocking, but that the effect is ten 
times stronger. 

2.5. The predicted mass evolution of clusters in the 
solar neighbourhood 

The decrease of mass due to the combined effects of stel- 
lar evolution, tidal stripping, spiral arm shocks and GMC 
encounters can then be described as 
dM _ (dM\ (dM\ (dM\ (dM\ 



Fig. 1. The mass evolution of a cluster with an initial 
mass of 10'* Mq in the solar neighbourhood. The mass loss 
due to the four separate effects is indicated. Encounters 
with GMCs are the dominant dissolution effect in the solar 
neighbourhood. 



1 1 



(9) 



tidal \ / sp \ / GMC 

with the terms given by Eqs.ElElEla'ndlHl We have solved 
this equation numerically for clusters of different masses. 
The results are shown in Fig.Qlfor a cluster with an initial 
mass of 10* Mq. The figure shows the total mass loss as 
well as the mass lost by each mechanism independently. 
Encounters with GMCs are the dominant dissolution ef- 
fect in the solar neighbourhood, contributing about as 
much as the three other effects combined. 

Figure 121 shows the ages of clusters when their remain- 
ing mass is and 100 Mq as a function of the initial mass. 
The almost linear part from log(Mi/MQ) = 3.5 to 5 has a 
slope of about 0.67. The figure also shows the dissolution 
times by the Galactic tidal field, predicted by BM03 for 
clusters with an initial concentration factor Wq = 5 in a 
circular orbit at i?o = 8.5 kpc. Our predicted timescales 
are about a factor 5 smaller, which agrees with the empir- 
ically determined dissolution time (L05). 

3. Comparison with observed age distribution of 
clusters in the solar neighbourhood 

Given the initial mass distribution of the clusters, their 
formation rate, CFR(t), and the time it takes for a dis- 
solving cluster to fade below the detection limit, we can 
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Fig. 2. The predicted dissolution times of clusters in the 
solar neighbourhood due to the combined effects of stel- 
lar evolution, tidal field, spiral arm shocks and encounters 
with GMCs, as a function of the initial mass. Full line: to- 
tal dissolution time. Dashed line: time when the remaining 
mass is 100 Mq. Squares and dotted line: dissolution time 
due to stellar evolution and the Galactic tidal field only, 
predicted by BM03. Cross with error bar: the value of ^4 
empirically derived by L05. 



predict the distribution of observable clusters as a func- 
tion of age or mass. L05 have derived an expression for the 
general case of a cluster sample that is set by a magnitude 
limit. Here we are interested in the prediction for a clus- 
ter sample that is complete down to a mass limit of 100 
Mq , because this is the ma ss limit of the unbiased sa mple 
of clusters within 600 pc of iKharchenko et al ] 1I2OO5I) (see 
L05). 
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For a constant CF R and a power law c luster IMF with 
a slope of —a = — 2 l|Lada fc Ladall2003() the number of 
clusters with M > 100 Mq as a function of age is 



(10) 



where Miim(i) is the initial mass of clusters that reach 
M(t) = 100 Mq at age t. (Clusters of age t with a smaller 
initial mass have M < 100 A'Iq by now.) Mmax is the max- 
imum initial mass of the clusters that are formed. The 
constant C is related to the star formation rate (SFR) in 
bound clusters as SFR = C ln(Minax A'^iim) for —a = —2. 

Figure 13 shows a comparison between the observed 
age distribution of clusters with M > 100 Mq within 
600 pc (from LOS) with the predicted distribution for 
Afjnax — 3 10** Af©. This value of Mmax is adopted because 
the observed distribution shows a steep drop at log t ~ 9.5 
(with only one cluster in the last bin) and Fig. [3 shows 
that this corresponds to Mi = 3 10* M©. The predicted 
relation for twice higher or lower values of M^ax agree 
worse with the observed relation. However, see discussion 
in Sect. 4. We have also calculated the expected age dis- 
tribution in case there was no mass-radius relation for the 
clusters, i.e. for A = 0. The downward slope of the result- 
ing distribution (not shown here) is significantly less steep 
than the one for A = 0.1 and does not fit the observed dis- 
tribution. 

The flattening of the predicted distribution at the low 
age end is due to the fact that clusters with an initial mass 
in the range of about 100 to 300 Mq quickly reach 100 
Mq (see Fig. 121). The bump in the observed distribution 
aro und log(t) 8.5 is due to a local starburst (see L05 
and IPiskunov et "all l)2006|) ). Notice the good agreement 
in the shapes of the predicted and observed distributions! 

The vertical shift that is applied to the predicted curve 
to match the observed one gives a value of C = 10^*'^^ in 
Eq. 221 which corresponds to a SFR of 4 10^ Mq Myr"^ 
for bound clusters in the range of 10^ < M^/M© < 3 10* 
within 600 pc from the sun. 

4. Discussion and summary 

We studied the dissolution of star clusters in the solar 
neighbourhood due to four effects: stellar evolution, tidal 
stripping, spiral arm shocks and encounters with GMCs. 
For this study we adopted the descriptions of GAPZ06 
and GPZB06 for the dissolution of star clusters by spiral 
arms and encounters with GMCs. We found that the last 
effect plays a dominant role in the solar neighbourhood. 

The cluster dissolution time due to spiral arms and 
GMCs depends on the density of the clusters, i.e. on 
M/r-n^ . This implies that the dissolution time is t^is ~ M^ 
with 7 = 1 — 3A if the radius of a cluster depend s on its 
mass a s ~ M^. We adopted A ~ 0.1 as found bv lLarsenI 
(Hooi) for clusters in spiral galaxies and so 7 ~ 0.7 for 
dissolution by both spiral arms and GMCs. This value is 
very similar to 7 = 0.62 predicted for dissolution by the 
tidal field only (BM03) and empirically derived for cluster 
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SFR = 4 10^ Mo/Myr 
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Fig. 3. The observed age distribution of an unbiased sam- 
ple of clusters with M > 100 Mq in the solar neighbour- 
hood within 600 pc (Karchenko et al. 2005; L05) in units 
of nr yr~^ is given by squares with the Poisson error bars. 
The full line shows the predicted distribution for a cluster 
sample with a maximum mass of A/max = 3 10* Mq and 
a SFR of 4 10^ Mq Myr^^. The dotted hue is for M^^^ = 
1.5 10* Af^ and the dashed line is for Af^,, ^ 6 10* Mo,. 



samples in four galaxies bv iBoutloukos fc LamersI l|2003|) . 

If there was no mass radius dependence for clusters in the 
solar neighbourhood, e.g. A = 0, then the predicted age 
distribution would have a shallower slope than shown in 
Fig. [31 since more old (massive) clusters would have sur- 
vived. 

Our calculated dissolution times of clusters in the so- 
lar neighbourhood are about a factor five smaller than 
predicted by BM03 for clusters in the tidal field of our 
Galaxy, with stellar evolution, binaries and two-body re- 
laxation taken into account. This is reminiscent of the 
short dissolution time of clusters in the central region of 
the interacting galaxy M51, where the empirical dissolu- 
tion time is even ten times shorter th an can be explained 
by stellar evolution and tidal fields l|Gieles et al. 2005ll . 
GMCs severely limit the lifetime of clusters in that galaxy 
also (see the discussion in GPZB06). 

The steep drop in the observed age distribution at 
f ~ 2 Gyr can be explained by an upper mass limit for the 
initial cluster mass in the solar neighbourhood of about 
3 10* Mq. However, this value is uncertain because it de- 
pends crucially on the completeness of the used sample at 
ages above 1 Gyr. (The sample contains only six clusters 
older than 1 Gyr.) The mass versus age distribution of our 
adopted sample, shown in Fig. 8 of L05, suggests that the 
lower mass limit of the observed clusters increases steeply 
for clusters older than 1 Gyr. Since the predicted value 
of dN/dt at any age depends on A/um and A/max as given 
in Eq. 1101 and Afiim increases when the lower mass limit 
increases, an increase in this limit implies an increase in 
Afmax derived from the observed age distribution. Based 
on this argument and the small number of clusters older 
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than about 1 Gyr in the observed sample, we conclude 
that the derived value of Mmax = 3 10* should be 
considered as a lower limit of the real maximum initial 
mass. 

The vertical shift applied to the predicted age distribu- 
tions to match the observed one indicates a star formation 
rate of 4 10^ Mq Myr'^ in bound clusters of Mi > 100 Mq 
within a distance of 600 pc, corresponding to a surface for- 
mation rate of 3.5 10~^'^ MQ-yr^^pc^^ . This is a factor 2 
to 3 smaller th an the SFR derived f rom the study of em- 
bedded stars bv lLada & Ladal ( 20031) because many of the 
stars are born in unbound clusters that dissolve within 10 
Myr. 

The very good agreement between the predicted and 
observed age distribution of clusters shows that dissolu- 
tion of clusters in the solar neighbourhood is dominated 
by encount ers with GMCs, as was already suggested by 
lOortI l)l958() . In fact, the good agreement may be slightly 
fortuitous because we have underestimated the dissolu- 
tion by two-body relaxation (see §1) and slightly overesti- 
mated the dissolution by encounters with GMCs, because 
we adopted the midplane density of GMCs whereas clus- 
ters may spend a faction of their lifetime above or under 
the galactic disk. Both effects are expected to be smaller 
than a factor two and may partially cancel out. 
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